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Abstract. We study the limiting behavior of the zeros of the Euler polynomi- 
als. When linearly scaled, they approach a definite curve in the complex plane 
related to the Szego curve which governs the behavior of the roots of the Taylor 
polynomials associated to the exponential function. Further, under a conformal 
transformation, the scaled zeros are uniformly distributed. 



1. Introduction 

Eighty years ago Szego |l8| studied the asymptotic behavior of the roots of the 
Maclurin polynomials associated with the exponential function. He found that if 
the roots are linearly scaled relative to the degree then the roots approach a curve 
S (see Figure 2) in the complex plane given by z G C such that = 1 and 

\z\ < 1. The behavior of the roots and poles of the Pade approximants and other 
Taylor polynomials have been analyzed ifTTIl . 

On the other hand, given any sequence of polynomials {p„(x)}, where Pn{x) is 
of degree n, asking how are the zeros of Pn{x) distributed in the complex x-plane is 
too general to get a reasonable answer. The best we can hope for is to focus on a 
special family of polynomials where a definite answer is possible. 

In this paper, we initiate the study of the asymptotic behavior of the roots of Euler 
polynomials E„ (x) which are defined by means of generating functions as 

(1-1) ^ = E£"«^- 

This generating function is listed among the principal generating functions by Louis 
Comtet Q for combinatorial applications. Their linearly scaled roots approach a 
curve related to the Szego curve S together with an interval on the real axis. 

Polynomials of binomial type were introduced by Rota and MuUin [7|. The rea- 
son for the name is that if 

n>0 



where D{x) is a polynomial, it follows that 
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which is a reminiscent of the binomial theorem. Herbert Wilf has the same de- 
scription in his book "generatingfunctionology" 1 12J . Strictly speaking, En{x) is 
not a polynomial of binomial type. However, it is close to being binomial type and 
has the simplest D{x) function. In this case £>(^) = ^. We hope to investigate the 
zero attractors of this wider class of polynomial families in the future. As evidence, 
the family of Bernoulli polynomials, for example, are easily handled with the tech- 
niques in this paper. A study of the behavior of their real zeros was recently done 

m. 

Let {pn{x)} be a sequence of polynomials. A set A in the .x;-plane is called the 
zero attractor of zeros of {pn{x)} if the following two conditions hold: 

a) Let Ag := Uyga ^{^^ £)' where B{x, £) is the open disc centered at x with radius 
e. That is, Ag is just the £-neighborhood of the set A. Then there exists an integer 
nQ{'£), for all n > hq, all zeros of Pn{x) are in Ag. 

b) For all ;c G A and for all £ > 0, there exists an integer ni(x,£) and a zero r of 
the polynomial p„j {x) such that r eB{x,z). 

Condition b) simply says that every point of A is an accumulation point of zeros of 

{Pn{x)}. 

The Euler polynomials En{x) are defined in dl.ll) . Since the nearest singularity to 
the origin of -^—^ are ^ = ±71/, it is easy to see that for all x G C, the power series in 
(11.11) converges absolutely and uniformly on any compact subset in |^| <n. In other 
words, although the polynomial En{x) is defined for all complex x but the power 
series is convergent only for ^ with |^| < n. By the Cauchy residue theorem, we 
have: 

n\ 27l//|^|=i (e^+l)^«+i ^ 
This integral expression is valid for all jc G C. Let x be replaced by nx and we can 
write the above equation as: 

The goal of this paper is to study the zero distribution of the polynomial E„(nx). 

2. A Generalization of the Szego Approximation 
We state two generalizations of the Szego approximation. Let 

Sniz) := t 



zJ 
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Proposition 1. Let S be a subset contained in \z\ > I so that the distance between S 
and the unit circumference |z| = 1 /s 5 > and let a be chosen so 1/3 < CL < 1/2. 
Then 

where the constant in the big O term is uniform for all z& S. 
Proof. By residue theory, we have for > 

S„-i{z) = — f y 

Note that ^ = z is a removable singularity of the integrand. Therefore, the above 
expression is valid for all complex z. For asymptotics in the region |z| > 1 + 5, we 
choose the contour to be the circle |^| =n. Thus 

Since nz is not included inside the contour |^| = n, a simple application of Cauchy's 
Theorem gives: 



:i J\r\=i l-z 



2Ki J\Q=l C 

where In ^ is the principal branch with —71 < arg ^ < 7t and n is a positive integer. 

Now we apply the saddle point method to construct the asymptotics of the integral. 
Since the critical point is ^ = 1 , the neighborhood of ^ = 1 must be carefully ana- 
lyzed. Let T| = with l/3<a< 1/2. The contour integral is decomposed into 
two integrals: 



1 r e«(C-lnQ 



2%iJ\Q=i ^-z 

where h is the integral on the circular arc in a small neighborhood of 1 : — T| < 
arg^ < r|, = 1, and the integral h is along the path in the remaining part of 
the circle. Set ^ = e'® in Ii. The Taylor expansion of the integrand in a small 
neighborhood of 6 = is worked out below: 
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Inserting these estimates in /i and carrying out some simplifications we get 



1 n /2(l + C>(ni-3a)) _ ^. 
— / 7; le civ 

e e ^^+^^^ ^^ (1+0(0)^6 



Since l/3<a<l/2, we have 3a — 1 < a. So the error term 0{n is absorbed 
into Hence, for |z| > 1 + 5, we see 

where the big O term holds uniformly for |z| > 1 + 5. If we put n0^/2 = m^, we get 

where (O = \J ^^^~2^ • Since a < 1/2, (O tends to with n. But, as n ^ oo, ^ = 

may therefore replace the limits of integration by ioo 
without altering the error term 1 + 0{n^^^'^). This gives 



\^27r(l -z) V n j-cx, y v27r/7(l-z) 

after taking /^^ e^" = into consideration. To justify that I\ gives the major 
contribution we obtain an upper estimate for I2 '. 

1 f e"5?(Q 

2% Jc 

where C is the contour determined by T] < |arg^| < tt, and |^| = 1. Obviously, 
9t(Q < cosTj. Note that |^ -z| > 5 for all zES. Hence 

Upon using 

costi = 1-^ + 0(ti4), 

we see that e«cosri ^ e«e"5"'"^"(l + But the factor g^i"'"^" = o(^). 

Consequently, 

h = o{h). 

This completes the proof of Proposition[l] . □ 

Next, the following proposition states the asymptotics of S„{z) in the region 
9t(z) < 1. 
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Proposition 2. For 1/3 < a < 1/2, we have 

e"^ v^(l-z)^ ^ 

where the big O constant holds uniformly for an arbitrary compact set ^ C 9^(2) < 
1. 

Thus the ordinary Szego approximation is generalized from the open disc |z| < 1 
to the open half plane 9^(z) < 1. 

Proof. The proof is actually very similar to that of Proposition [T] Here we use a 
suitable integral representation for (z) in the region. We start off with: 

'"-'^'^ = 2^ifcZ^z^'^^ 
where C is any closed contour encircling the origin. This integral representation is 
valid for all complex z. We insert nz for z to obtain 

1 /■ - {nzT 

2-KiJcC,-nz C," 

Let K be an arbitrary compact set in the domain 9^(z) < 1. We choose the circular 
integration contour C : C, = —R + {R+ —n < 9 < TC with R so large that K 
becomes strictly included inside C. Let z belong to K and replace the contour C by 
nC. This gives 

After a change of variables ^ — > and the use of the Cauchy Integral Formula, we 



2ni Jc^-z{nQ'^' 



This implies 

e"^ 27rz Jc ^ — z 

where C is the circle ^ = —7? + (i? + l)e'^, — tt < 9 < tt. The critical point is still 
^ = 1 which is a point in the contour C. It is clear that the saddle point method can 
be applied to the contour integral. The procedure is very similar to what we did 
in the previous case. We omit the details. This is how we prove the statement of 
Proposition|2l □ 

We can rederive the same result by using the integral representation of in 
the original derivation of the Szego approximation. Szego [9J used 

, n'-'e- r 
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to obtain the result of the proposition in the disc |z| < 1. A good elaboration of his 
approach can indeed offer a proof of PropositionEl The saddle point method as we 
presented in the proof is just an easier way to get the result. 

Proposition 3. (Jet Wimp) (Uniform Approximation ofS„{nt)for t >0) 

Sn{nt) 



uniformly for t >Q, where ^(?) 

5(0 

2 

andforallt,Erfc{t) := e^'^ ds. 



1 1 /2 

\t —I — In? I ' , and 

1, for 0<?<1, 
0, for t > 1. 



This version of the Szego approximation was proved by Jet Wimp (personal com- 
munication). 

3. A Decomposition for the Euler Polynomial 
The Euler polynomial En{nx)/n\ is decomposed as the sum of two polynomials: 
Proposition 4. Let /j be an integer > and let F^il,) be given as 



1 



Then we have 



■I 

k=0 



1 



{2k+i)Tii{^-{2k+i)Tir 

EJnx] 



+ ■ 



1 



-{2k+l)Tli{^,+ {2k+l)Tli) 



n\ 



Mn,^,{x)+Kn,/x) 



where 



H,^. (t) 



2E 

^=0 



Sn-\{nx{2k+ l)ni) ^ S,j-i{—nx{2k+ 1)7^0 



((2/t+l)7r/)«+i (-(2/t+l)7l/)"+i 
Here 5„(z) := Y!j=oZ-^ / jl the n-th partial sum ofe^ as usual. 
Proof. The following integral representation for E„{x) is valid for all x G C: 



En{x) 2 f 



n\ 271/ J|^|=i (e^ + l)^«+i 
Let X be replaced by nx to get 



En{nx) 2 r 



,4' 



d^ 



n\ 271/ j^^|=i I ^ / ^(e^+i; 
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Note that for each integer k>0, 

1 



{2k + l)ni{^-{2k+ 1)71/) {-{2k + l)ni){^+{2k+ l)m) 

is the sum of singular parts of ^(^J_^_i'^ {2k+ l)ni and —{2k+ 1)71/. Hence Ffj{^,) is 

analytic in the annulus < |^| < (2/i + 3)71. The bigger // is, the bigger the domain 
of analyticity of Fj_i{b,) is. Hence 



(3.1) x|£ 

<k=0 



1 1 



{2k+l)m{^-{2k+l)Tii) {-{2k+l)m){^+{2k+l)Tii 
A typical term of the above sum is: 

2 f e^^\ 1 



2niJ\i^\=i \ % I {2k+l)Tii{^-{2k+l)%i) 



1 r -i^y 



271/ (2A:+ l)7l/7j^|=i (2A:+ 1)71/ ;'^ (2fc+ 1)71/ 

after using a geometric series expansion. Obviously, the series is uniformly conver- 
gent on IQ = 1, we carry out the integration term by term. By the Cauchy integral 
theorem only those terms with I <n — \ survive. Thus we get 

1 (-1) «- 1 / 1 ^' 



7i/((2A:+l)7t/)2^'^ \{2k+l)'Kij Jj^| 

71/ {{2k+ l)7i/)2 ;to^ {2k+ 1)71/^ {n-l-l)\ ^ ' 
- {{2k+l)%i)^to (n-l-iy. ^^^"^^^'"'^ 



(-2)((2.+ l)7./)-«-£M^^±^ 
j=o J- 
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where 5„(z) is the n-th partial sum of e^; that is, = E^^qzV jl. In a similar 
way, we obtain 

2 r f e'^^\ " 1 
2^4l=iV W {-{2k+l)ni){^ + {2k+l)Tii) 

Inserting these back into (13.11) . we complete the proof of the proposition. □ 

Proposition m is important since it provides asymptotics for the Euler polynomials 
in various regions. The asymptotics for Mn,^(x) (for any fixed /j) can be easily found 
by the classical saddle point method. It is 

where the big O term holds uniformly for £< \ l/x\ < (l/j + l)n + e. The unifor- 
mity can be justified using the fact that F^{^^) is analytic in < |^| < {2ij + 3)n. 
The asymptotics for Kn^^{x) can be obtained from the generalized Szego's approxi- 
mation. Using these asymptotic approximations we will prove that the point set K 
defined below is the zero attractor of the Euler polynomials and can also determine 
their density distribution. 

4. The Zero Attractor 

Let the point set K be defined by the graph in Figure 1 : 
The point set K consists of the curves as indicated and the real interval [— 
K is symmetric with respect to the imaginary axis in the .jc-plane. We will show 
as the first step that K contains all accumulation points of the zeros of En{nx). In 
the second step we will carry out the density calculations of the zeros residing in 
an immediate neighborhood of K. As a consequence, this will establish that every 
point K is an accumulation point of zeros of En{nx). Hence K is precisely the zero 
attractor of E„{nx). The density calculation gives also the statistical distribution of 
zeros which includes the information about the fraction of zeros along each segment 
of K. The work of proving that K contains all accumulation points is divided into 
the three parts in the following lemma: 

Lemma 1. (a) Let \xq\ > ^, then xq is not an accumulation point of zeros of the 
Euler polynomials En{nx). 

(b) There is no accumulation point of zeros in the region \x\ < except for real 
numbers. 

(c) Ifx is a non-real accumulation point of zeros in the region 3^ < < then we 
must have either 

Uji/e^-^' 1 = 1 or I -xTiie^^"^' 1 = 1. 
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Figure 1. 



Proof. Let us prove part (a). Now 



where /(^) = ^ — In^. We choose the principal branch of In^ here. To invoke the 
saddle point method, we need to find the critical points which are roots of /'(^) = 
or 1 — ^ = 0, that is, ^ = 1 is the only critical point in question. Observe that 

|exp/(^) I = ^ = e^°^^, where ^ = e'^. Now e^°^^ attains its maximum at = 0, 
hence at ^ = 1 . By the saddle point method |I3 we have 

En{nx) / / \ 1/2 




\/27m giA+l 

where the big O constant can be made uniform for a given compact set K C {x : 
\x\ > ^}. From this, we see that for any given set K C {x : \x\ > ^}, there exists no 
such that for all n>nQ the polynomial En{nx) has no zeros in K. Hence the proof 
of part (a). 
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Next we prove part (b). Assume that there is a non-real accumulation point of 
zeros in the region \x\ < 3^, say xq. Of course, xq 7^ 0, then there exists an integer 
/Jo > 2 and a positive number £ > such that 



(2^/0+1)71 



+ £ < l^ol < 



1 



(2^/0-1)71 



Since xq is an accumulation point of zeros of E„{nx), there exists an infinite se- 
quence of integers nj such that 

Enj {n jXnj) = and jc„ . xq as j -> 00. 

We may assume that for all large j, {xnj} is in the region 



+ £< \xo\ < 



1 



+ e, 



{2ijo + 1)k - (2^/0-1)71 

so we apply Proposition|l|with /j chosen as /jq — 1 and keep /jq fixed in the following 
arguments. Use the asymptotics of Mn^^{x) we see that En ■ {njXn ) = implies 



(4.1) 



2 jxnjep 

n ./JTjXnj 



1 



1 



where the big O term holds uniformly for all x in the region 

1 1 



(2^/0+1)71 



+ £ < be < 



+ £. 



Note that the summation in Kn^^Q-i{xnj) is running from to //q — 1 and for all large 

A typical term in (jc„ .) is 

Snj-iinjX„j{2k+ 1)71/) 

Introducing t = Xnj{2k+ l)ni. For < k < /jq — 2, we have \t\ < I —e. Hence we 
apply Szego approximation of Proposition |2| Thus 



e"Jt 



Snj-l{njt) _ 

{{2k+l)%ip+^ ~ {{2k+l)mp+^ 

'-7ss=T^<""''"'('^''<"; 



l-3aN 
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The big O constant in the above approximation is uniform since \t\ < ^^jpy + £• 
Introduce the function 

g(x) := or = (Tixie^-"^')-^. 

enxi 

Thus the above equation when expressed in term of g{x) becomes 
Snj-ijnjt) _ iexn^rjg'yii2k+l)xnj) 
{{2k+l)Tii)"j+^ ~ i2k+l)ni 



(4.2) X 



where l/3<a< 1/2. Since xq is not real, so we may assume that xq is in the lower 
half plane. Thus xq = roe^'^" with < Go < and (2^,)^ 1)^1 + £ < ro < j2jJ^T)n- 
Now, 

gTlrosinOo g-irrosineo 



l^(^0)| = > \g{-XQ) 



enro enrQ 

^7trsin6Q ^ . ^ . . . . l 

Furthermore, ^— — , as a function of r attams its minimum at r = „ a ■ For each 
1 < k < ijQ — 2 (remember /jq > 2) we have 

ro<{2k+l)ro<{2iuo-l)ro<-< ^ 



K 71 sin 00 
Hence 

(4.3) 1^(^0)1 > \g{{2k+\)xo)\ > \gi-{2k+\)xo)\. 
Also, since xq lies in the interior of the Szego's domain 

|2oe^"'"| < 1 and |zo| < 1- 

Here, zo = Tlxoi. So, 

l^(^o)|>l. 

This prepares the behavior of the terms in K„^ijQ-i{xnj) with <k <fio — 2. The term 

5„ . _ 1 (mjX„ . (2jUo — 1 ) m) 

corresponding to ^ = /jo — 1 in the sum Kn,/jQ-i{Xnj) is ■' 1^ .^„^._|_i — which 

can be estimated using Proposition |3l Let = y^^.x;„^ (2/Jo — 1)^^ and apply the 
proposition to get 

\Snj-i{njXnj{2iJo-l)-Ki)\ = |5„ ._i ((n 1 )?^J | 

<5„._i((nj-l)|?„,|) 

(4.4) < e(-^-i)k^ol 2./| ^^|^^"P Erfc {./J^l^iM)) (l+0{^)) . 

Note that \t^^^\ < ( ^,^^^'_^-^^ +e)(2^o - 1)7T = 1 + (2„) - l)7r£ as nj oo. Recall that 

g{xo) = e^'^'^'/enxoi, 
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where 

xq = Ixole''^", < 00 < 7t 

and 

1 1,1 

-j — — 7T-+e< N < Tw' A'o > 2. 

(2//o + l)7t (2^0 -1)7^ 

A detailed study of the Szego's curve, defined as = 1 and |z| < 1, shows the 

following features (see Figure 2): 




ego curve 



Figure 2. 

point B is the intersection point of the curve with the imaginary axis and equals 
i/ e, point C is the intersection point of the curve with the negative real axis and is 
—0.278 • • • (the unique real root of e^~^ = —x) point A is the intersection point 
of curve with the circle |z| = 1 /3. 9t(A) ^ —9.861 x 10~^ that comes from solving 
the real root of 1/9 — x^ = e^^^^ — x^ for x. Hence, zo = tixqI falls in the interior 
of the Szego's curve (also to the right of the imaginary axis in the z-plane). This 
implies that 



-zo 



Recall ^(|?^|) = 



1 — InUuQ 1 1 So it is easy to see that 



0(1) 





^MO 


) 


ho\ 








1 
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and Erfc( a/ — 1^( | ) ) := / e ^ ds = 0(1) uniformly. Hence from 

(1^ we get 



(4.5) 



for some absolute constant K. Now we are ready to see a contradiction from (14. It . 
Dividing (14.11) by {x„je)"j we get 

y!^'^;^^ (i + o(l)) =0 



mated as 



Snj-i{njXnj{2iAo - 



{eXnji2iuo~l)TliYJ 



(4.7) 

Similarly, 
(4.8) 



< K 



K 



Snj-l{-n jXnj ( 2/Jo - 1 ) TTZ ) 



(-e.^„.(2A/o-l)7r/)"> 



'nj) 'Kn,lJ{)- 




- 1 in {exnj) 








e |x„ . 1 (2^0 










1)71 j 






\xnj 1 (2;U0 - 


1)71 j 



By (14.21) and (14.31) the dominant term in {ex„j) "■'K„^ijq-i{x„j) corresponding to sum- 
mation indices < ^ < /jq — 2 is the term corresponding to = 0; that is, the term 



Tti \ ^/2nnj 1 -XnjTli 
The other term with the same index = is 
W^^A 1 1 



-.g-"J{xnX'^ + 0{n'-n) 



-HI \ ^ /Innj 1 + XnjTii 
which is still dominated by the above term. Note that 



:^-">(-x„.)(l+0(nl-3a)) 



TtX„ ■ I 

e J 



enxnji 



ttLt,, sine„ 

en \x„i I 



where Xnj = "■' ■ Since Xnj —^xo = rQe^'^" with < Gq < 7r, this implies 

sin0„^. > for all large nj. Hence, \g{xnj) \ > — r — r implies 



eTt\x„ 



(4.9) 



1 



enlxn: 
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Note also that for /jq >2 implies ^ > j2~iyK- comparing (14.91 ) with (14.71) 
and (14.81) we see the dominant term corresponding to ^ = still dominates the term 
corresponding tok = /jq—I. Consequently, we infer from equation (14.61) that the left 
hand side becomes arbitrarily large as nj oo. Hence the left hand side is non-zero, 
contradicting to the right side of the equation. This completes the proof of part (b). 

We now proceed to prove part (c). We show that if xq is a non-real accumulation 
point in the region < |jc| < ^, then either 



\xo%ie 



1 — X()7t( I 



1 or 



-XQllie 



l+XQlt/ I 



That ^ and real x such that ^ < |x| < ^ are points of the attractor is a consequence 

of the density calculation. Since by assumption xq lies in < \x\ < i, we can 

certainly choose £ > sufficiently small such that .x;o is in + e < < ^ — e. 
Again as in the previous cases we may assume 



Xq = \Xo 



e"''^",whereO<0o<7r. 



Also the same £ works for an infinite sequence of zeros Xnj in 

£ such that Xn - ^ Xq as nj ■ 



1 , , 1 

— -^£< JC < - 
371 



7C 



Now use Proposition m with the choice = 0. So, 



where 



Mnflix) 



KnAx) 




Sn-\{—nxni) 



and 
(4.10) 



■ + 



+ ■ 



1 



^(e^ + 1) 7r/(^-7r/) {-%i){^ + m 



Now Enj(njX„j) = implies M„.o(x„.) +KnjQ{xnj) = 0. Since Fo(^) is analytic in 
the region £ < |^| < 371 — £, the asymptotics for A/„^.y(x„.) is obtained as: 



(4.11) M^.^Xn.) 



[Xn/PPQ 



njXfij 



where the big O term holds uniformly. Note that since Xnj xq we have ( +£)7r < 
|j:„^.71/| < — £)7i:. So XrijTii lies in a compact set in the half plane 9tx < 1. Hence, 
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we can invoke Proposition|2lfor the asymptotics for Knj q {xnj). Combining this with 
(I4.11I) in (I4.10I) and expressing the result in terms of the function g{x), we obtain 

71 X„. y^Xnj V "7 

ni V ^ v^(l -Xnjjni V J 



as Hj 00. 



Note that there are terms in y ^Fo{\ /x„j)-^^j^ which are to be cancelled in the 

above equation and the order term 0{^) is absorbed in This observation 

leads to a simplification: 

(4.12) 

Also note that |^(^o)| > !,?(— -^o)| (strictly greater). This implies 

Therefore, if |^(xo)| > 1, then the left hand side of (14.121) becomes arbitrarily large 
in modulus for large rij. This is a contradiction. But if |g(xo)| < 1, then the left 

hand side of (14.121) goes to iAo_|_| ' ^ oo, a non-zero number which is still 
a contradiction. Hence we must have 

k(^o)| = l 

that is, 

|;co7r/ei--'^°'''| = 1. 
This also shows that zq '■= x^^i is a point on the Szego curve: 

|zo| < 1, Izo^'^^^^l = 1, -7^/2 < argzo < 7^/2. 

that is, 9?zo > or zo lies in the shaded segment of the curve or Xi^ is in the rotated 
Szego's curve. 

To show there is no real accumulation point xq satisfyin g ^ > > we can 
still use the asymptotics that leads to (14.121) . In particular, (14.121) still holds. Note 
that in this case Xn^T^i lies in the exterior region of the rotated Szego curve in the 
.)c-plane. This implies \g{xrx^ I < 1- Hence in the limit as rij ^ oo in (14.121 ) we get 

^ 1 _^ 
7tglAo-(-l ' 

a contradiction. This finishes the proof of part (c). □ 
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This establishes the fact that the point set K contains all points of the zero attrac- 
tor. 

5. Determination of the density of zeros 

The next step is to carry out the density calculation. The Euler polynomials 
satisfy 

£„(i-x) = (-i)"+i£„w 

We have 

E,{nx) = (-1)"£„(1 -nx) = (-l)"£„(-n(x- 1/n)). 

From this it is easy to see that the zero attractor does has the reflection symmetry 
with respect to the j-axis. The Euler polynomials are polynomials with real coeffi- 
cients. It obviously has the symmetry with respect to reflection about the jc-axis. To 
describe the density we choose the lower half of K for discussion. 




Figure 3. 

The image of the points A, 5, and C in the .x;-plane under the mapping of z = nix 
are denoted by A' ,B' , and C' respectively. Further, the subsequent images A" , B", 
and C" are the images of A', B', and C' under C, = ze^^^ (see Figure 3). 

We show: 

Theorem 1. The image of the zeros ofE„(nx) along the arc BC in the x-plane are 
uniformly distributed in the C,-plane along the corresponding circular arc B"C". As 
a consequence, the fraction of zeros residing in a neighborhood of the arc BC in the 
x-plane is 

Tl/2-\/e _\ 1 
27t ~ 4 2Ke' 

Theorem 2. The real zeros ofE„{nx) falling in the line segment AC are uniformly 
distributed in the segment AC. As a consequence, the fraction of zeros residing in 
the segment AC is ^. 

Before enumerating the zeros, we mention some well-known facts about the 
Szego curve: The mapping C, = ze^^^ is conformal L4J in |z| < 1. The Szego curve 
is defined by: 
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l^gl-^l =1, < 1 

The mapping ^ = ze^^^ is in a neighborhood of z = 1. It is 2-to-l in a neigh- 
borhood of z = 1 since ^ = and ^ 7^ 0. Therefore, to enumerate the 

number of zero images inside a contour we will face a difficulty of inverting the 
function in a neighborhood of 1 in the ^-plane. For this reason we will choose 
a contour that does not enclose 1. But how many zero images are left out in a 
neighborhood of 1? We will show that it is of order o{n). To prove this we use 
Jensen's inequality. However, the use of Jensen's inequality requires a knowledge 
of the function En{nx) in a neighborhood of —i/n. Fortunately, this knowledge is 
sufficiently provided by Proposition[T]and Proposition |2| 

Theorem 3. (Jensen's Inequality) Let h{z) be analytic in the disc \z~a\ <R and 
let < r < R and m be the number of zeros ofh(z) in the disc \z — a\ <r. Then we 
have the following inequality: 

R\"' ^maxi,_a\=R\hiz)\ 



r J 



A proof can be found in most books of complex analysis. Let denote the 
number of zeros of En{nx) which lie in the disc l^t^ — ^ | < e in the x-plane. 

Proposition 5. For all sufficiently small £ > 0, there exists no(£) such that for all 
n>nowe have 

Yi"^<j^('^ln(l + 0(£))+ln| 
where the big O constant and K are absolute. 
Proof. Since does not vanish in j-'c — ^ | < £, the function 

(5.1) hn[x) .= — 

n\{ex)" 

is well-defined and has the same number of zeros as En{nx) in j.'c — i| < £. We 
apply Theorem|3to hn{x) on the disc — ^| <2e. Thus 

/2£y^"' . max|^_i^|^2El^«W 



2Ye = n < 



We need asymptotic estimates for | /?„ ( ^ ) | and max | ^_ i, | ^2e I (^) I • Fi^^t, for | /?„ ( ^ ] 
we apply PropositionHlwith /j = 0. Thus 

En{nx) 



n\ 
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where 



and 



The asymptotics of M„ o(^) is known: 
1 



5„_ 1 {nXKi) ^ Sn- 1 ( —nxni) 



i-ni) 



n+l 



(5.2) 
Recall 



nfi 



Til 



2 / / e \ « ^ , ..TO 
Fo{%i 



Kn,0 — 



n \ \Tii 
2 



Sn-l{n) + 



l+0(- 



(-1) 



ii+l 



We will show that ^„,o(^) is the dominant terni for En{n^^. Use Proposition 13 to 
get (with t = n/{n — \) and n ^ n — 1 ) : 



(5.3) Sn-i{n)=e" -2_^Evfc{V^K) j (^l+^^y^ 
where^=|;^-l-ln^|^/2 = -15^(1+0(1)). So 



Recall that 



Erfc(x) 



e ^ ds 



e ^ ds— I e ^ ds 







+ 0{x) 



Erfc 



as X ^ 0. Hence 

Erfc(v^^r^^) 



Inserting these estimates back into (15.31) we have 



V 2 Vn - 1 n 
1 , 



2 



5„_i(n) 



l+0( 



1 



Similarly using Proposition 13 we have 

Sn-i{-n)=0{e-"). 
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These estimates give 



^«,o(-|:) = (-).)(^rfi+o(^; 

Til ni ni V \/n 



Comparing (15.21) with the above we see that the order of ^n,o(^) is small than that 
of Knfii^) by a factor of y^. Hence 



n\ Ki Til 

and 

We now estimate maxi \hn(x) I . To this end we use PropositionHJwith ^ = 0: 

En{nx) 



where 




The asymptotics of Mnfi{x) is: 



Sn-i{nxTli) Sn-ij-n) 

(71/)"+! (-71/)"+! 



V TT \ X y'nx / \ n 

Care must be exercise to handle the asymptotics of Sn-\{nxTli) . Recall that x is on 
the circumference — I =2e. 

I Til I 

Choose two points p and q (see Figure 4) such that p is the mid-point of the 
circular arc between the circle |;c| = ;^ and the horizontal tangent line to the circle 
\x\ = ^at the point x = The point q is similarly selected. Note that as e — > 0+, 
the distance between p and the circle \x\ = ^ is of order 0{e^). This is so because 
the horizontal line y = ^ is tangent to the circle \x\ = ^. The image of the arc 
qtp under the map z = xTli is a circular arc in the z-plane where Proposition 1 is 
applicable. Thus for all x on the arc qtp we have 

(5.4) Sn-x{nxTli) = (i +o,(;,(i-3a) 

\'2Tln{xTli - 1) ^ 
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Figure 4. 



where stands for the £-dependence for the big O constant. Taking (15.41) into 
consideration we get 



hu (x) 



Vo(-)-(l+0(i)) 
Tl X X n 



+ 



+ 



2g"(x) 
ni 



2%{x%i- \ ] 



2>g"(-^) 
{-ni) 



-'\-x) 



-(l+0e(nl-3a)) 



where the third term is obtained from an analogous application of Proposition El to 
Sn-\{—nxni). We now estimate hn{x) as follows. The first term 

^Fo(i)i(l+O(i) 
n X X \ n 

is obviously < K, an absolute constant. The second term 



2g'\x) g-''{x) 



2n{xni — 1 



< 



K 



;i+0e(nl-3a)) 



IxTli — 1 1 



K 



< -(l+C>,(ni-3«)) 
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The third term 



1 - 



-(l+0e(nl-3«)) 



\/2nn{l +xni) 

because —xni lies outside of the Szego curve and \g{—x) | < 1. Hence for x on the 
arc ptq on — ^ | = 2e we have 

|/^.WI<|(l+0e(nl-3a)), 

where K is an absolute constant. When x in the arc psq on the circumference | ■'c — ^ | = 
2£, we invoke Proposition|2l Thus 

2V^g"{x) 



K{x) 



X X 



^o(-)-)(l+0(-)) + 



X 



1- 



-(l+0e(nl-3a)) 



V^nni 1 — xni 



(5.5) 



i-ni) 



1 



^/2nn(l +xni) 

The magnitude of g{x) for x on p^'f^r is estimated as follows. Recall xni = 1 + lenie'^. 
This implies 

,-2£jre'(''+'t/2) 



1 — xJl; I 



Ixni 



Ixn 



This yields 



1 



xnie 



Therefore, we infer from (15.51) that 



\hn{x)\<K 



< 



1 -2£7r 



,2e7i 



1 -2e7r^ 

Note that in deriving the above equation the fact that < 1 was still used. 

Combine these estimates to get 

„2e3t 



max \hn{x) \ < K 



U-i =2e 



1 -2£7l 



Recall that 
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yi"^ < In 



This implies that 

^) +0,(„>/»«))J/,„2. 

This finishes the proof of Proposition |5l □ 

We now come to determine the density of zeros. First of all, we refer to the 
mapping relation in Figure 3. In general, we let N„{a, (3) be the number of image 
points of zeros of En{nx) in the ^-plane that fall in the angular sector a < arg ^ < p. 
Now let an arbitrary be given in the interval 

^ ^ ^ 1 

O<0< 

2 e 

Our goal is to establish the following proposition: 

1 
Proposition 6. lim -A^(O,0) = — . 

Proof. Recall that A^„(0, 0) is the number of image points in the ^ -plane that fall in 
the angular sector < arg^ < 0. We alleviate the problem that the straight edges 
of the above sector may contain some image points by perturbation. The reason for 
requiring that no image points of zeros fall on the straight edges is to guarantee an 
application of the argument principle. Since the Euler polynomials are polynomials 
of rational coefficients, the roots are algebraic numbers in the x-plane. This implies 
that the totality of image points in the ^-plane is countable. Hence we can choose 
an arbitrary small number £ > so that the straight edges of 3£ = arg ^ avoids all 
images points of zeros. The following arguments assume a fixed £ > and n will 
be ^ oo eventually. It is now obvious that we have the following inequality: 

(5.6) A^„(O,0) -iV„(3£,0 + £;j = O (y^")) , 

where a null sequence £', > is chosen so that no image points are on the straight 
edge arg ^ = + £^. Define 

_ En{nx{Q) 

^"^^^ •= ;^iMQF^- 

Note that: (1) the image points are zeros of fn{Q and (2) fn{Q = hn{x{Q), where 
hn{x) was defined in (15.11) . Here the function x{Q is the inverse map of the map 
from the x-plane to the ^-plane. x{Q is 1-1 except in a small neighborhood of 

We introduce the contour e^g/^ defined as (see Figure 5): Fi U r2 U U r4: 

re'3^1-£/2<r< l+£/2, 
(l+£/2)e'^3£<(^<0 + £;„ 
^e'(e+e;,),l_e/2<r< l+£/2, 
(l-£/2y'^,3£<^<0 + £;. 




EULER POLYNOMIALS 



23 




Figure 5. 



We apply the argument principle to enumerate the number of image points. Thus, 

-dC, 



A^„(3£,0 + <) 



J_ r t^MQ 



(5.7) 



-(3(/i)+3(/2)+3(/3)+3(/4)) 



where /, = /p., 1 < ?' < 4. We shall handle h first. If we pull the contour r2 back in 
the z-plane, it falls in a compact set in the half plane 9^(z) < 1 so that we can use 
Propsition|2lfor the asymptotics of 5„-i(«z). Thus with an application of Proposi- 
tion |4| (with /J = 0) we have 



En{nx)^ 



+ 



+ 



Recall 



n XX 
ni 

i-ni) 
1 



l+0{-) 
n 



v27rn( 1 — x%i 



-{l + 0,{n'-^^)) 



v27m(l +x'Ki 
1 
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SO 

1 X 1 1 

°^x'^~eiA + l %i{l~xTii) {-ni){\+x%i)' 
Canceling some common terms to simplify the expression we get 

En{nx)^/n _ [2 1 
n\{exY "VTTeiA+l 

(5.8) I 1 + 0(i) + + S'(--)^"(f"+n ^ o(„,-3„; 

1 n Ki 

Furthermore since the pulled back x's lie outside the zero attractor so \g{—x) \ < 
\g{x) \ < 1. In this way we see that 

£,(nx)v^eiA+l ^ , 

lim — ; — > 1, uniiormly 

n\{ex)" a/S/tt 

or 



hn{x) ^177^' uniformly 



or equivalently 



hn{x{Q) ^^(Q— , uniformly for ^ G 

& I J- 



By a theorem of uniform convergence of analytic functions we get 

hMQ) ^ eiA(0 + i 

uniformly for ^ G Recall f„{Q = hn{x{Q). Hence by integrating the above 
along r2 we get 



1 



This implies 



1 1 X .1 

--3(/2) = 0(- 



Let us handle U now. The asymptotics in (15.81) still hold good in this case. But since 
X now is inside the zero attractor so 

\g{x)\ > 1 and still \g{x)\ > \g{-x)\. 
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Rewrite (15.81) in the form 

,l-3a 



)) 



(5.9) +v|Z^-,^(^,(_,))«. 

Note that ^ = zei"^ and |^| = 1 -e/2 since ^ G Also, \t,g{-x)\ = \g{-x)/g{x)\ < 
1. Hence we see that 

uniformly and we get similarly 

uniformly for ^ G Integration along r4 gives 

Note that the orientation of r4 gives the correction of the "— " sign. Hence, 

3(/4) =n3 (/(0 + £;,-3e)) +o(l). 



This implies that 



Caution must be exercised when we handle /i ( can be analogously handled). 
Although ( 15.91 ) still holds good for /i, but |^| varies from 1 - e/2 to 1 + £/2. It 
is not clear how ( 15.91) becomes useful to us. We now employ a number-theoretic 
argument and Theorem|3to obtain a useful estimate for /i . The argument we follow 
is inspired by that of K. Chandrasekharan |IT]. Let / be the number of points on 
Fi so that 9t(/„(Q) = 0. We insert these points into Fi and decompose /i as the 
sum of integrals whose end points are two consecutive points where 9t(/„(Q) = 
except possibly the beginning integral and the final integral. A typical integral to 

consider is of the form 3(/j' y^J^) , where a and b are two consecutive roots of 

5K(/„(Q) = along Fi . By a change of variable, we have 



UQ 7 \Jc, i 

where the contour C„ is the image of the segment ab under the map C, — > fn{Q- We 
comment that C„ intersects the imaginary axis at its two endpoints and at no other 
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interior point of C„. No matter what, we always have by Cauchy's Theorem 

J C„ ^ •/ Kai, ^ 



where Kab denotes the semi-circle with segment f,j{a)fn{b) as the base so that Kab 
lies in the same half plane as C„. Note that /^^^^ db,/^, is either 0,in, or —in. Hence 



K-ab ^ 



<7r. 



It follows that 
(5.10) 



|3(/i)|<(/ + l)7r. 



Here we have + as an upper estimate because of the possible inclusion of the 
beginning and final integrals. To estimate Z, we shall use the Jensen's inequality: 
Define 

hniQ ■= 2 ' 

where hn{x) was defined in (15.11) . Let us note the following two properties of /z„(Q: 

(1) hniQ is an analytic function of C,. 

(2) When C, is real and 1 -£/2 < ^ < 1 +£/2, then ^e'^^ G Ti and x(^e'^^) and 
further — jc(^e^'^^) are complex conjugates in the x-plane. 

Hence 

The reason we choose 3£ instead of £ is for convenience as the sequel will show. 
Now we regard h„{Q as an analytic function defined on the disc |^ — (l+£/2)|<2£ 
in the ^-plane. (This is so because the circle {C,e^^^ : |^ — (1 +£/2)| = 2£} does 
not include ^ = 1 in its interior, for otherwise, x{C,e'^^) would not be well-defined 
in the {^e'^^ : |^-(l-f£/2)| < 2£}.) Note that the disc : |^-(l+£/2)| < £} 
contains the real interval [1 — £/2, 1 +£/2]. Each root of 9^(/„(Q) along Fi is a 
real root of hn{Q in [1 — £/2, 1 -|-£/2]. Let / be the number of roots of /j„(Q in 
{C,: 1^ — (1 +£/2) I < £} (possible complex roots are counted in /), then obviously 

(5.11) 1<T. 

Apply TheoremEin the disc : |^ - (1 + £/2) | < 2£} to get 



(5.12) 



2£ 
£ 



2'< 



max|^_(i+e/2)|=2e 


hniQ 




/l„(l+£/2) 





Recall the definition of /i,, ( 1 + £/2) : 



hn{l +£/2) = - [h„{x{{l+e/2)e''')) +h„{-x{{l +e/2)e~''% 
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The point C,= ( 1 + e/2) e^^^ lies in a region in the ^-plane so that the pull back z( ( 1 + 
e/2)e^^^) lies in the half plane 9t(z) < 1. Therefore, we can use (15.51) to determine 
the asymptotics of /j„(x((1 +e/2)e^'^)). Similarly for hn{-x{{l + £/2)e^^'^)). In 
this case we have 



^(-x((l+£/).-3'B))| ^ U(^((i+e/2),^-))| < 1 



3ie\ 



SO 



(5.13) 
where 



lim /z„(l+e/2) 
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A study of the curves defined by9t( ^,y| ^ ) = 0, that is, e^/-^ + 1 = 5 G M leads to 

1 



X 



In V\+s^ + / arg(— 1 + si) 

These curves represent points in the x-plane where e^/^ + 1 = purely imaginary 
numbers .We have branches of the curve which is a consequence of the multi-valueness 
of arg( — 1 +si). These curves cluster at;c = 0. The point =^((1 +£/2)e^'^) is in 
a small vicinity ofx= —i/n. Note that 

dx —i 



ds 



2 • 



This means that the curve has a vertical tangent atx= —i/n. In the ^-plane 



d 
de 



[\+£/2)e 



3ie 



+ 3z. 



e=0 



It is easy to see that the pull back of Xe does not lie on any of the branches of the 
curve 

1 



91 







Moreover, since e""' + 1 = 0, it is not hard to see that 



1 



K 

^7' 



for some positive constant K. This implies from (15.131) . 



K{l+e/2) 



K 
> —. 

£ 



Next we estimate 



max 

|C-(l+e/2)|=2e 



max 

|C-(l+e/2)|=2e 



En{x{Q)V^ 



n\{ex{Qy 
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we can still use (15.51) for the purpose. This is so because the puUback z{Q is still in 
< 1. Proceeding similarly and using ( 15.51) , we get 



max 

K-(l+e/2)|=2e 



h„iQ <i^[(l + 0(£)rv^ + 0e(nl-3«)] 



where the big oh constant in 0(e) is an absolute constant. Further, combining (15.111) 
and (15.121) . we get 

/<j^ln[(l + 0(£)rv^ + 0e(ni-3«)]. 

Insert the above in (15.101) , we have 

K 
\n2 

Similarly we obtain 



3(/i)|<— ln[(l + 0(£)rv^+0e(ni-3«)]. 



Now from (15.61) 



|3(/3) I < In [(1 + 0{e)rV^+0,{n'-'^)] 



-A^„(O,0) - V„(3£,0 + e;j = -0(4"^: 



n n n 

Inserting the estimate from 3 (/,) , 1 < / < 4, in the above and taking the limit sup, 
we get 

hmsup ^ ' ' < — + (9(e). 

fi — >oo yi 27r 

An analogous lower bound can be obtained, that is, 

. ^ A^„(O,0) ^ 
liminf — ^ > 0{e). 

Since e can be made arbitrarily small, we conclude 

lim 1a^,(O,0) = ;^. 

This ends the proof of Theorem [T] □ 

We next calculate the density of zeros in the interval [0, ^] . The strategy will be 
the same to that of Theorem[T] However, the technical details are slightly different. 
We outline the steps below. 

Lemma 2. (a) For every £ > 0, let be the number of zeros ofEninx) in the disc 
{x : \x\ < £}. We show that < ^{Ke{1 + 0{e))"). 

(b) For allQ <b < we construct the rectangular contour Ci U C2 U C3 U C4 = C 
as shown in Figure 6. 

Let Nn{£,b) be the number of zeros ofE„{nx) enclosed inside the contour C. In 
general, we introduce the notation: let N„{a,b) be the number of zeros of En{nx) 
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Figure 6. 

enclosed in the rectangular contour with vertices a — id, b — /5, b + /5, a + id. Then 
for alley and for all sufficiently small 5 > we have 

5, 1 
£' n\n2 



-N„{E,b) <b + 0(5) + 0(£) + 0(tan-i -) + -l-ln(i^e (1 + 0{£)f). 



A consequence of the above two lemmas is: lim„^oo j^N{0,b) = b. Thus the real 
roots of En{nx) are uniformly distributed in [— - 



_i_ j_ 

Tie ' Jle J 



Proof. Let us show part (a) first. Since ^"^^"^-^ is analytic on the disc {x : \x\ < 2£}, 
we apply Theorem |3] to get 

(5.14) 



(I) 



2^^ < 



max|^ 


=2e 


E„{nx) 




£«(0) 













En(O) 



can be obtained from the integral representation: 



En{0) 



1 



1 



n\ 7rz7|^|=i (e^+l)^"+i 

Here caution must be exercised because for even n except n = 0, E„{0) = 0. This 
means that E„{x) has a zero at x = for all even n>2. Therefore, we must modify 

-E (nx) 

the function if we still want to use Jensen's inequality. We consider , " when n 
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is even. In this case, an equivalent Jensen's inequality is: 



2N. < 



maxi 



x|=2e 



Enjnx) 



limx^o 



-E„{nx) 



We will show the former inequality has the main features in the proof, while the 
latter can be treated similarly. From integral representation (n is odd now) we apply 
the Darboux method [5|. In this case the nearest singularities are ±ni. One shows 
that 



En{0) 



1 



n 



(-1) ^ (-1) 



1 .t: . 1 ^ 



(5.15) 



(7r/)"+i ' {-Tii)"+^ 
The corresponding quantity when n is even is 



+ o(— ) 



lim,,^o ]^En{nx) _ 1 ^.^ 1 / 2 / 



Now we know that f|^|=i (ei;+i)^"+i = (n is even here). We use the above to 
rewrite limx^o -£^n('^-''^)- Thus 



lim -Enijix) 

x^Q X 



lim i) 



n\ x" Ini Jj^|=i (e^ + l)^"+i 

' 2 i 



n\ 2niJ\^\=\ l)^«^ 
2 / 



^ i) — 7 — dt {n even) 

(n-l)!2m/|^|=i (e^ + l)^« ^ ^ ' 



which is still of the same feature as in (15.151) . In this way we see it does not matter 
whether n is odd or even, they can be handled in a similar way. Now we consider 



max 

U|=2e 



En{nx) 



We use Proposition |4| with /Ji chosen sufficiently large so that ^2iii-^\)n < 2£ < 
(2,ui-i)jt " Here the additional assumption that the arbitrarily small number It is 
not of the form ,^ iiw ?^ ^ 0, will not hurt the arguments we present here. Then 



with /J = /Ji — 1 in PropositionHlwe get 

En{nx) 



-l{x) +Kn^^^-l{x). 
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The asymptotics for M„^^j_i (x) is 

71 



X ^nx J \ n 



and the asymptotics for comes from applying Proposition |2| This is so 

because \x{2k+ \ )Tli\ < \x\ (l/ji — 1)71 for all < < /Ji — 1. When \x\ = 2e, we 
find \x{2k+ < 2e{2/ji — 1)tz < 1. This means the Proposition |2l is applicable. 
Furthermore, the largest term comes from k = 0. Thus 



En{nx) 



(5.16) 



= (x.)"u/|iV,i-i(-)^fl + 0(- 
V 71 X ^nx \ n 

+ ^^(l+o(l)) + ^f^(l+o(l))] 



Note that on = 2£, x is in the interior of the rotated Szego curve ^xnie^^-"^' | = 1- 
Hence \g{x)\ > 1, -<:)| > 1. On the upper semicircle of |.)c| =2e, \g{—x)\ > \g{x)\, 
while on the lower semicircle of \x\ = 2e, \g{x) \ > \g{—x)\, and on the real axis 
\g{x) \ = \g{—x)\. No matter what we always have 



e \x\ n 



l+6>(e) 
e \x\ n 



,for \x\ = 2e. 



Similarly |^(-jc)| < Hence from (l5l6ll we conclude 



E„{nx) 



< 



K^{\ + 0{z)Y 



for some absolute constant depending only on £. Using (15.141) and (15.151) . we get 

< 



max|;r 


=2e 


En{nx) 




£„(0) 





Taking logarithms we get 



u\+o{z)r 



N,<—\n{U\ + 0{z)Y. 



This ends the proof of part (a). 
We now prove part (b). 

N{z,b) 



271 \Jc hn{x) 



dx 
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where h„{x) = ^"}"(lly^ ■ Decomposing C into Ci,C2,C3, and C4 as shown in Figure 
6, we get 

(5.17) A^(e, Z;) = i-3 (7i ) + ^3 {J2) + ^3 (73) + ^3 (74) , 

where 

= [ ^Jjc, l</<4. 



Let us focus on 7i first. We infer that (15.161) still works for all x on C\ provided 
Z{x) =5 is sufficiently small. This is so, for \g{x)\ = = ^^^j^^, where for 
X eCi,x= \x\ < < 7r. This observation leads to 

that is, g'\x) is the dominant term among the terms in Kn,fn-i{x). Hence h„{x) 
g^"{x) — * ^ uniformly. This implies 



hn{x) X 
Integrating the above along Ci we get 



h' (x) 1 
" -n{iii ) -^0, uniformly 



ni{b-io-e + io) - In 



/5 



1 ^ „ , b-id 
Tl[b — £) — arg 



that is. 



271 V "e-/5 
^ + 0(5) + 0(tan-i?) 



lim -^3 (7i) = J + 0(5) + 0(£) + 0(tan-i -) 



Similarly, we obtain 

\ b S 

lim —7(73) = - + 0(5) + 0(E) + 0(tan-i -) 
n^oo 2%n 2 £ 

(In this case, g"(—x) becomes dominant, rather than g"{x) .) An estimate for ^3 (72) 
comes from the observation that the change of arguments for g{x) on the vertical 
segment b — id, b is of order 0(5) (15.181) . Hence, 
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Similar to (15 .141) we can likewise prove that 



2nn n\n2 
Inserting all these estimates into ( I5.17I) we get 

1 5 1 

-Nn{e,b) <b + 0(5) + 0(e) + 0(tan-i -) + -— -ln(^e (1 + 0(e))"). 
n E R In 2 

Again taking (15.141 ) into account we get 

-Nn{0,b) < Z7 + 0(5) + 0(E) +0(tan-i -) + ^ln{K,{l + 0{e)r). 
n E n In z 

This implies 

1 5 

limsup -A^„(0, b)<b + 0(5) + 0(e) + 0(tan-i -) + In (1 + 0(e)) . 

We get a similar lower bound for 

1 5 

lim inf -A^„(0,Z7) > - 0(5) - 0(e) - 0(tan"i -)- In (1 + 0(e)) . 

n^'x n E 

But the above is true for all E > and all 5 > 0, hence 

lim -Nn{0,b) = b. 

In this way we have established all claims we have made in this paper. □ 
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